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Abstract. Results on dissipative isoscalar modes of a hot and dilute nuclear droplet are presented. As 
compared to the adiabatic limit (part I), realistic dissipation yields a substantial reduction of the growth 
rates for all unstable modes, while the area of spinodal instability in the (£>,T)-plane remains unchanged. 
The qualitative features of multifragmentation through spinodal decomposition as obtained in the adiabatic 
limit are not significantly affected by dissipation. 



PACS. 21.60.Ev Collective models - 21.65.+f Nuclear matter 
tation - 25.70.Pq Multifragment emission and correlations 



25.70.Mn Projectile and target fragmen- 



1 Introduction and summary 

In the preceding part I Q| we have introduced a collective 
model which allows to study the eigenmodes of a nuclear 
droplet as function of its density g and temperature T . 
We refer to this paper for a thorough discussion of the 
model and related publications. The description is based 
on the diabatic approach to dissipative collective motion 
and - in the local density approximation - yields equa- 
tions of motion for small amplitudes, where the mass and 
stiffness tensors are obtained analytically. The model is 
suited to explore systematically characteristic properties 
of hot nuclear droplets as function of their densities. While 
this is interesting by itself and particularly so in the region 
of the liquid-gas phase transition, such studies also yield 
basic information on multifragmentation by spinodal de- 
composition, and hence may give insights which are com- 
plementary to dynamical simulations. 

In Q we have studied isoscalar modes in the adia- 
batic limit, which is defined by a vanishing relaxation time 
(r = 0, instantaneous intrinsic equilibration). However, 
the adiabatic (or thermodynamic) limit is quite unrealis- 
tic for nuclear systems, because reasonable values for the 
relaxation time are of the same order as the characteristic 
times of collective motion in the region of the liquid-gas 
phase transition. Therefore, we extend our study here to 
arbitrary values of the relaxation time, and thereby in- 
clude dissipation. 



By varying r between infinity and zero our diabatic 
approach allows us to treat the continuous transition be- 
tween the two elastic limits, i.e. the diabatic limit(r — > 
oo), where the dynamical distortions of the local Fermi 
sphere are not destroyed by two-body collisions, and the 
adiabatic limit (t = 0), where the Fermi sphere is assumed 
to be restored instantaneously. 

The method of solving the secular equation for the 
eigenmodes, as well as the determination of the relaxation 
time r and the evaluation of the diabatic part of the stiff- 
ness coefficients (not included in Q) are presented in sec- 
tion H When expressed in terms of the temperature, the 
relaxation time of a Fermi-liquid is practically indepen- 
dent of the density g, yielding rates H/t ~ T 2 /3MeV. In 
addition to the adiabatic spinodal we find also diabatic 
spinodals for the bulk modes of infinite nuclear matter 
and nuclear droplets, however not for surface modes. 

Since the secular equation for the eigenmode energies 
is of third order, there exist three roots for < r < oo 
instead of two for the adiabatic and diabatic limits. As 
functions of r the eigenvalues move around in the com- 
plex w-plane on generic trajectories which depend on the 
positions in the (g, T)-plane (cf. section 3). 

For realistic values of the relaxation time r the stable 
and unstable regions in the (g, T)-plane are examined in 
section |] and discussed in relation to the adiabatic (ther- 
modynamic) limit of [Q. The main effects of dissipation 
are summerized as follows. 
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— The spinodal lines for bulk and surface modes are not 
changed by dissipation. 
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— However, the growth rates of unstable bulk and surface 
modes are quite sensitive to dissipation and, due to the 
sensitivity of the relaxation time r on temperature T, 
are reduced from their adiabatic values by factors 1/2 
to 1/4 in the region of interest. 

— The qualitative differences between soft and stiff equa- 
tions of state (EOS), with larger spinodal regions and 
larger growth rates for the stiff EOS, survive, when 
dissipation is included. 

— Dissipation leads to some increase in the splitting of 
growth rates for modes with different multipolarities / 
and number n of radial nodes. With decreasing density 
the quadrupole mode I, n = 2, becomes unstable first, 
followed in sequence by the higher multipoles I — 3, 4, 5 
and nodal number n = 1,2. 

The results show that the qualitative features of multifrag- 
mentation by spinodal decomposition survive when dissi- 
pation is included (section 5). 

In order to keep the presentation concise we have to 
refer frequently to sections, figures and equations of part 
I jj] and denote these references by a prefix I. 



2 Determination of dissipative eigenmodes 

As discussed in jl| , the eigenvalue equation for the collec- 
tive modes in harmonic approximation reads (eq. (1.19)) 



B xx > + C' xx , 



lu + i/r 



Cxw = 



(1) 



with the mass tensor B = {B\\>}, the adiabatic stiffness 
tensor C = {C\\> } and 



= {ci A ,} = c-; 



the difference between the diabatic (C) and adiabatic stiff- 
ness tensors. The relaxation time t characterizes the decay 
time of deformations of the local Fermi sphere. 
By substituting u> = ifl one can rewrite (|l|) as 

/2 3 + -i7 2 +B" 1 (C + C')^+-B" 1 C = 0, (3) 
r r 

where all coefficients are real numbers. This problem is 
equivalent (see Q for details) to the diagonalization of 
the nonsymmetric real matrix 



I 



- x C/r l-^C + C) I/r 



(4) 



All results presented in section 4 are based on this diago- 
nalization procedure. The mass tensor B and the adiabatic 
stiffness tensor C are given in jl| . In the following subsec- 
tions we present expressions for the relaxation time r and 
the difference C of the diabatic stiffness tensor with re- 
spect to the adiabatic one for compressional and surface 
modes. 



2.1 The relaxation time r 

According to |3) the relaxation time r in the nuclear Fermi 
gas with equal numbers of neutrons and protons is given 

by 



3.25 avp(g) ge* 



(5) 



where sp, a, vp , g and e* denote the Fermi energy, nucleon- 
nucleon cross-section without Pauli blocking, Fermi veloc- 
ity, nucleon density and the excitation energy per particle, 
respectively. The factor 3.25 results from numerical eval- 
uations of collision integrals. Since we are dealing with 
small amplitude vibrations, the energy stored in the de- 
formation of the thermal Fermi distribution is negligible 
compared to the thermal energy, i.e. 



' £th 



7T 2 T 2 

4e F (g) 



, , g eq \ 2/3 m*(g) 



a T' S2t 
' Q 



m*{g eq ) 



(G) 



where £p{g) oc g 2 ^ 3 m*(g eq )/m*(g) and the reduced level- 
density parameter a$ = clq/A for stable nuclei with equi- 
librium density g eq has been introduced. Replacing e* in 
(||) by the expression (^|) we find 



'(8) 



a T 2 m*(g eq ) 



with the p-independent and dimensionless constant 



(3tt 2 /2) 1 /3 



1 



6.5 



2/3 
<7 geq 



= 0.36 



(7) 



(8) 



(2) 1 



where a = 40 mb (4 fm 2 ) and g eq — 0.14 fm~ 3 has been 
used. Further realistic values for 77 have been reported in 
and range from 0.25 to 0.52 for ao — 0.1 McV -1 . Note 
that the range is roughly [0.15. .. 0.30] if the Fermi-gas 
value 0.067 for ao is applied [jj. Through m* the relax- 
ation time depends only weakly on the density g. In the 
numerical calculations of stability and instability (section 
4) 77 = 0.25 is used, hence we have for the relaxation rate 
typically H/t rs T 2 /3MeV in the spinodal region, where 
m* (g) / m* (g eq ) rj 1.2. 



2.2 Stiffness tensor C 

According to (||) C denotes the difference between the 
diabatic and adiabatic stiffness tensors. Since their inter- 
action parts are identical, C ' reduces to the difference 



= C (T) - C 



{T) 



(9) 



of the diabatic and adiabatic contributions from the in- 
trinsic kinetic energy. In evaluating this difference in the 
local density approximation we start from the correspond- 
ing expression for the diabatic intrinsic kinetic energy 



, fi 2 V'<^(r')V'<ft k (r') 
2m*(f?(r')) 



(10) 
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Fig. 1. Adiabatic (solid lines) and diabatic (dashed lines) spin- 
odals in the (g, T)-plane for infinite nuclear matter (indicated 
by oo), gold-like (A=197, Z=79, JV=118) and tin-like (A=W0, 
Z=N=50) nuclear droplets calculated for a soft EOS (SkM*). 
Here, go = 0.16 fm -3 is the normal nuclear matter density. 
Note that the equilibrium density of the droplets is around 
0.8520- 

for protons and neutrons, separately. Here r' = r + s(r) 
is the position shifted by the displacement field s(r) and 
^k(r'), q(t') denote, respectively, the distorted normal- 
ized plane- wave functions and density at the point r'. The 
momentum distribution /(k) remains fixed in the diabatic 
limit and is normalized to the density g of the unperturbed 
homogenous droplet. With (I. A. 3) and (I. A. 8) expression 
(|l0|) in the local-density approximation reduces to 



zpintr 
^kin 



lr<R 2m* (g) 



(11) 



where g is given by (LA. 12). The distorted momenta k are 
determined according to (I. A. 3) by 



dr'. 



with r'j = Tj 



s 3 (r), and hence 
(k 2 ) = J ^ m ^ = \^{p-. 



(12) 



(13) 



where the isotropy of f(k) has been used. Inserting the 
expansion 



dr' n 



= 6, 



E 

k 



dsi 



dr. 



(14) 



up to second order in s, we obtain with (1.41) 

(k 2 ) dwb = (k 2 ){l-|div S + J2(did 5 w) 2 } (15) 



1,3 



in the diabatic limit. 

With the adiabatic expression for (k 2 ) from @ 

(k 2 ) adlab = (k 2 ){l - hivs+ |(divs) 2 + 1 Y^ 9id l w ^ 



hi 



(16) 



(note the misprint in (I.C.2)) we finally can evaluate C 
according to (H) as 



C' xx , = e d 3 r 



h 2 (k 2 ) 

2m*(g) 
d 2 1 2 



(17) 



dqxdqx 



E(^ w ) 2 -g( divs ) 



where the density dependence of m* does not contribute, 
because the bracket is already of second order in q\. 



2.2.1 Compression modes 

The displacement fields for the compressional modes have 
been denned by (1.34,1.41) and discussed in detail in 0], 

s(r,i) = Vw(r,t) = qx(t)V X x(v) . (18) 

A 

with the normalized functions 

Xx={nlm}(r) =Afnlji(K nl r)yj n (n) , (19) 

where ji(K n ir) denote the spherical Bessel functions with 
ji(nniR) = 0, and yj 71 ^) the real spherical harmonics. 

By applying the relations (I.C.4) t o ( I.C.7), noting a 
wrong sign in (I.C.6), we obtain from (|l7j) the final result 



A A' 



i— n,p 



2(2 
3 1 3 



'V 



R 2 



Sii'S n 



(20) 



where the contributions for protons (i =p) and neutrons 
(i =n) are summed. 

(q-) 

When comparing the diagonal elements C' xx with C xx , 
i.e. with the contribution from the intrinsic kinetic en- 
ergy to the adiabatic stiffness coefficient (cf. eq. (1.46)), 
one realizes that C' xx is of almost the same magnitude, 
also always positive, because (K n iR) 2 3> 1, and hence 
a major contribution to the diabatic stiffness coefficient 
Caa = C' xx + C X \ , cf. figs. 1.3 and 1.4. As a consequence 
the diabatic spinodal in the (g, T)-plane, which essentially 
is determined by C X \> (g, T) = with I — 2, n = 0, lies well 
inside the adiabatic spinodals as shown in fig. [f]for infinite 
nuclear matter and two nuclear droplets with A=197 and 
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100. Like for the adiabatic spinodals the large finite size 
effect is due to the finiteness of the minimum K n i -value 
for the droplets. In the diabatic limit this almost leads 
to disappearance of the spinodal region. Furthermore, the 
diabatic spinodals for the gold-like and tin-like droplets 
are almost identical like in the adiabatic limit. 



2.2.2 Pure surface modes 

In accordance with j5) and (1.59) we define the displace- 
ment field for pure surface modes by 

S M)= E n^V(r<3T), (21) 

l>2,m 

which is consistent with incompressibility of the flow up to 
first order in the collective variables Qi m . This is sufficient 
for the evaluation of (|l7|) up to second order in Qi m , and 
hence with V • s = we have 

C' xx , = i4\q,T)\ J d'rj^mw f . (22) 

i—n,p i,j 

Taking the additional factor \J (21 + 1)/Att in the defi- 
nition of our collective variables into account, we find from 

§ 

/W-tM* (23, 

and 

2e^ (g ,T) ^- 1 ) 3 f + 1 ^ 3. (24) 

i—n,p 

For completeness we also give the adiabatic stiffness tensor 
for the surface modes 

C = C (S) + C (c) (25) 

with 

C (S) =C^^1(^) 2/3 , (26) 

\ Qeq / 

where C^f^' 1 denote the expressions (1.63) and (1.64) from 
|jTJ and g eq the density of the cold stable droplet. The 
g- and T-dependence have not been reported in Q. No 
diabatic spinodal exists for the surface modes because 

a > \c\. 



3 Generic features of dissipative modes 

For simplicity we consider a single mode such that the 
matrix equation (1) reduces to 

G^tuj) = -Boj 2 + C — ^—+C = 0, (28) 

where G(u>) denotes the Green function. Since the cou- 
pling of different modes does not alter the qualitative de- 
pendence of the three roots (poles of the Green function) 
on the parameters B,C',C and r, we are able to reveal 
generic features of the dissipative modes from (p8[). 

3.1 Motion of poles as function of r 

For different values of the mass and stiffness parameters 
we illustrate these properties by the motion of the roots 
(poles) in the complex w-plane as function of r, which 
is freely varied from the diabatic limit (r — > oo) to the 
adiabatic limit (r — > 0), thereby creating trajectories for 
the three poles in the complex cj-plane. For the presen- 
tation in fig. U we have chosen the compressional mode 
l,n = 2,0 of a gold-like droplet for the soft EOS (SkM*) 
at five points in the (g, T)-plane (cf. fig. |l|), ranging from 
the region of stable nuclei via the vicinity of the adiabatic 
spinodal around T = 3 MeV into the regions of adiabatic 
and diabatic spinodal instabilities. 

— The top part (a) of fig. ^ corresponds to a gold nucleus 
at small temperatures (C, C > 0). For r — > oo there 
exist essentially two poles of diabatic vibration. In this 
limit the third pole corresponds to an amplitude con- 
stant in time added to the time-dependent vibration. 
For increasing dissipation, i.e. decreasing r- values, the 
diabatic poles become complex, moving into the lower 
half of the complex cj-plane, and hence correspond to 
damped oscillations. The third pole (sometimes called 
thermal pole) moves also down on the imaginary axis. 
For t — » the two vibrational poles approach their 
adiabatic positions on the real axis while the thermal 
pole quickly escapes to — ioo. 

— When approaching the adiabatic spinodal (C > 0, C > 
0) around T = 3 MeV, we observe an interesting de- 
formation of the pole trajectories as illustrated in part 
(b) and (c). The vibrational pole trajectories approach 
each other and in (b) almost touch. Only somewhat 
closer to the spinodal line the trajectories are no longer 
separated, and hence the poles loose their original phys- 
ical meaning. Starting from the diabatic limit (r — > oo) 
the two vibrational poles meet on the negative imagi- 
nary axis and then move up and down (towards — ioo) 
while the thermal pole moves down slowly on the imag- 
inary axis. The upward moving "diabatic" pole collides 
with the downward moving "thermal" pole and finally 
approach the adiabatic pole positions as r — * 0. 

— The trajectories for a position somewhat inside the 
adiabatic spinodal (C < 0), but outside the diabatic 
spinodal (C = C + C > 0) are shown in part (d). Now 
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Fig. 2. Generic behavior of the eigenvalues of eq. (^) in the 
complex oj-plane as functions of the relaxation time r, pre- 
sented in (a) to (e) for different regions of stability and insta- 
bility (different values of C and C + C , see text). Limits for 
r —* and t — » oo are indicated. The bullets on the trajectories 
correspond to realistic r-values obtained for r\ — 0.25. 



adiabatic vibrational position in the lower half plane 
and — ioo, respectively. 

For a position inside the diabatic spinodal (part (e)) 
both C and C = C + C are negative, such that both 
pairs of adiabatic and diabatic vibrational poles lie on 
the imaginary axis. With decreasing r-values the insta- 
bility pole moves again upwards on the positive imag- 
inary axis, however starting now from the finite value 
in the diabatic limit t — > oo and reaching the adiabatic 
value for r — > 0. For decreasing r the thermal pole col- 
lides with the upward moving lower diabatic pole and 
after a detour in the lower complex half plane both 
poles collide again on the imaginary axis and finally 
move to the adiabatic position in the lower half plane 
and to —ioo, respectively. 



3.2 The pole of instability 

We look for the pole on the imaginary axis lo = i"f for 
finite values of r. According to (03) we have 



B-i 1 +C- 



7 



+ C = 



(29) 



7 + 1/t 

which for the dimcnsionless quantities a — C/C , [3 = 
l^B/C and 6 = Ty/C'/B reads 



a 



-f) 2 







P+1/0 



(30) 



This universal function is illustrated in fig. ||for (3 > —1/9, 
i.e. above the singularity (3 — —1/9. All curves enter the 
region of instability [3 > through the origin, where the 
adiabatic stiffness coefficient C, and hence a becomes neg- 
ative. The family of curves map the whole region between 




for all finite values of r there is a pole on the posi- 
tive imaginary axis reflecting the spinodal instability 
encountered in the spinodal region. The two diabatic 
poles move down with decreasing r and, after collid- 
ing on the negative imaginary axis, they approach the 



Fig. 3. Universal curves for the evolution of the instability 
poles /3 = jy/B/C as functions of a = C/C' for different 
relaxation times r between the adiabatic (r = 0) and diabatic 
(r -> oo) limits: 6 = t^/TFJB = 0, 0.3, 1.5, 6.3, 25.5, 102, 410, 
1638, oo. 
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Fig. 4. Growth rates (Im(u) > 0) of the unstable compres- 
sional quadrupole modes (7 = 2) as functions of density at 
T — 3 MeV for a gold-like droplet and the soft EOS. Since the 
coupling between the radial modes is small, the eigenmodes 
are still well characterized by the nodal numbers n = 0, 1, 2, 3. 



Fig. 5. Growth rates (Im(u) > 0) of the unstable compression 
modes with multipolarities I = 2, 3, 4, 5 as functions of density 
at T = 3 MeV for a gold-like droplet and the soft EOS. Note 
that the curve I = 2 is identical with the curve n = of fig. 



the adiabatic limit r = and the diabatic limit r — > oo 
(0 — > oo). In the unstable region between — 1 < a < the 
growth rates are quite sensitive to the relaxation time and 
vanish in the diabatic limit. 

4 Effects of dissipation on stability and 
instability 

Eigenvalues are obtained by numerical diagonalization of 
the matrix (^j) for given values of density and temperature. 
As pointed out in modes with different angular mo- 
menta I are decoupled. Coupled arc modes with the same 
I and different n (n is the number of nodes in the spherical 
Bessel function determining the radial density profile) . As 
the energy of modes grows substantially with n, while the 
couplings between the modes is relatively small, it is suffi- 
cient to take into account only a few modes with the small- 
est n-values. We have studied all modes with I = 2, 3, 4, 5 
and included in the diagonalization the radial modes with 
n = 0, 1, 2, 3 for each i-value. As mentioned in section 2.1 
we use 77 = 0.25 for calculating the relaxation time r. 

4.1 Compression (bulk) modes 

The results on the compression modes are presented in 
figs. I to |. 

4.1.1 Dependence on n and I 

Like in the adiabatic limit (cf. figs. 1.5 to 1.7) more and 
more modes with n > for a given Z-value become un- 



stable (lm(u>) > 0) with decreasing density as illustrated 
in fig. H for I = 2. However, dissipation reduces the mag- 
nitude of the growth rates by factors 1/2 to 1/4 and, in 
contradistinction to the adiabatic limit, no crossings of 
different modes as function of g is observed. Furthermore, 
the curves are not that well degenerate at small densities 
as in the diabatic limit, such that the modes with n — 0, 1 
dominate. This behavior is typical for all multipole modes 
as well as for the soft and stiff EOS. 

Figure H (the analog to fig. 1.6) shows the largest growth 
rates (Im(u>) > 0) for different multipoles (I = 2,3,4,5). 
With decreasing density the growth rates for different l- 
values come close to each other. However, again the de- 
generacy is not as perfect as in the adiabatic limit. Still, 
different multipolarities will strongly compete in the decay 
of the droplet. 

4.1.2 Dependence on g and T 

In contradistinction to 0] we present in the regions of sta- 
bility instead of the roots u> only the damping rates 7 = 
Im(cj) < of the vibrations. These quantities are relevant 
for the restoration of the equilibrium shapes. In the spin- 
odal regions (the regions of instability) the positive imag- 
inary roots determine the growth rates 7 = Im(cj) > of 
the unstable modes. The following discussion of stability 
and instability is based on these quantities. 

Figure H (the analog to fig. 1.7) shows the largest growth 
rates in the unstable (shaded) region and the smallest 
damping rates in the stable region for a gold-like droplet 
described by the soft EOS. The adiabatic spinodal line is 
not affected by dissipation, because it is essentially deter- 
mined by Uadiab for I, u — 2, 0, the coupling from finite 
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Fig. 6. Contour plot of largest growth rates in the unstable 
(shaded) region and smallest vibrational damping rates in the 
stable region for dissipative compression modes of a gold-like 
droplet described by the soft EOS. The numbers on the lines 
denote the rates Im(hui) in MeV. All modes with I — 2, 3, 4 and 
n — 0, 1, 2, 3 are taken into account. The dashed line indicates 
the diabatic spinodal. 

r-values between different n being smafl in the relevant 
regions g 0.25g eq , T ss 4 MeV. There the growth rates 
are reduced by factors 1/2 to 1/4 as mentioned before. 
Furthermore, a pronounced effect from the diabatic insta- 
bility is seen for small T (large r-values), consistent with 
the behavior of the growth rates around a — — 1 in fig. 0. 
Since r oc 1/T 2 the growth rates and damping rates ap- 
proach zero on the real axis to the right of the diabatic 
spinodal. 




1 i i i i i i i i i i I i i i i i i i 
0.0 0.2 0.4 0.6 0.8 1.0 
P/Po 



Fig. 7. Contour plot of largest growth rates in the unstable 
(shaded) region and smallest damping rates in the stable region 
for dissipative surface modes of a gold-like droplet described by 
the soft EOS. All modes with I = 2, 3, 4 are taken into account. 

small densities and temperatures, while for large tempera- 
tures and densities outside the bulk spinodal only surface 
modes are unstable. 



4.3.1 Dependence on size 

Figure ^ (the analog to fig. 1. 10) shows the results in 
the (g, T)-plane for the tin-like droplet. Like in the adi- 
abatic limit the compression instabilities of a finite nu- 
clear droplet depend only weakly on the mass. However, 
the surface instabilities are pushed to considerably smaller 
densities for the lighter system. 



4.2 Surface modes 

Figure (the analog to fig. 1.8) summarizes the results 
for the surface modes. As for the bulk modes, dissipation 
reduces the growth rates of instabilities by factors 1/2 to 
1/4 in the region of interest (T « 5 MeV). The spinodal 
line is again not affected by dissipation. The shift in fig. ^ 
to smaller densities as compared to fig. 1.8 is not an effect 
from dissipation but due to the inclusion of the correct 
density dependence also in the Coulomb part (25) of the 
stiffness coefficient, which was neglected in M. 



4.3 Combined plot of bulk and surface modes 

A combined plot of compression (bulk) and surface modes 
is presented in fig. || (the analog to fig. 1.9). Like in the 
adiabatic limit the compressional instabilities dominate at 



4.3.2 Dependence on EOS 

Figure [l0| (the analog to fig. 1. 11) displays the combined 
instabilities for the gold-like droplet and the stiff EOS. 
Dissipation does not alter the qualitative dependence on 
the EOS. The region and degree of compressional insta- 
bility is considerably larger for the stiff EOS. Again the 
growth rates in the relevant region are reduced by factors 
1/2 to 1/4 due to dissipation. As far as the calculations are 
comparable these results are consistent with those quoted 
in @. 



4.3.3 Sensitivity to dissipation 

Figure |ll| illustrates the sensitivity of the growth rates on 
the reduction of the relaxation time to half of its value. 
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Fig. 8. Combined bulk and surface instabilities (soft EOS) for 
a gold- like droplet. The dashed lines indicate the adiabatic and 
diabatic spinodals. 



Fig. 10. The same as in fig. fel but for the stiff EOS. 



> 

CD 




0.0 0.2 0.4 0.6 0.8 1.0 

P/Po 

Fig. 9. The same as in fig. H but for a tin-like droplet. 



While the surface instabilities remain almost unchanged, 
the bulk instabilities gain roughly 30% in the growth rates. 
Since a factor of 2 defines roughly the uncertainty range 
for the relaxation time we can consider the growth rates 
presented in figs. ^ to 10 reliable to about ±30%. 
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Fig. 11. The same as in fig. ^, but for smaller relaxation times 
with rj = 0.12. 



5 Effects on multifragmentation by spinodal 
decomposition 

In section 1.5 we have shown how spinodal decomposi- 
tion through adiabatic surface and bulk instabilities may 
explain qualitatively and to some extent even quantita- 
tively characteristic features of multifragmentation reac- 
tions. In discussing the effects of realistic relaxation by 
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two-body collisions we consider the following three con- 
secutive stages of multifragmentation: 

- the initial formation of a hot nuclear droplet with den- 
sity QfSS Q e g, 

- expansion into the spinodal region and 

- spinodal decomposion. 

The various reactions (cf. which are studied with 
respect to multifragmentation, yield nuclear droplets typ- 
ically with excitation energies E* /A w 6 MeV. At these 
high excitations two-body collisions lead to local intrinsic 
equilibration (relaxation of local Fermi sphere) on short 
time scales of r w 10 fm/c, cf. section 2.1. Thus thermal 
intrinsic equilibrium is attained already at the beginning 
of the expansion. However, this equilibration does not in- 
clude the relaxation of deformations which are produced 
in the formation of the heated droplet. As shown in figs. || 
and bulk and surface modes relax on the time scales 
of h/{2 MeV) « 100 fm/c and h/(l MeV) w 200 fm/c, 
respectively, and hence survive to a large extent the ex- 
pansion process into the spinodal region (duration time « 
(30 . . . 50) fm/c, cf. ||). Therefore, in order to study for 
example the bulk instabilities in the spinodal region, it is 
important to create the hot droplet in a gentle way with- 
out too large density distortions. Otherwise the surviving 
initial density fluctuations could mask the characteristics 
of spinodal decomposition. 

The expansion of the heated droplet is fast, as men- 
tioned before. However, no significant deviations from in- 
trinsic equilibrium is expected. The reason is that the ex- 
pansion is well approximated by the radial flow v(r) = ar, 
which keeps the density homogeneous within the droplet. 
This flow scales the single particle energies like the tem- 
perature with g 2 / 3 /m*(g) such that the product of single- 
particle energy and temperature, and hence the single- 
particle occupation probabilities remain constant. Thus 
starting from an equilibrium distribution, the system stays 
during this isentropic expansion in equilibrium even with- 
out relaxation by two-body collisions. Consequently, within 
the collective model described by (|l|), the diabatic and 
adiabatic stiffness coefficients are identical, and hence no 
damping of the expansion occurs (C = 0). Dynamical 
calculations || have shown that for not too large initial 
excitations (E* /A < 10 MeV) the expansion reaches turn- 
ing points at temperatures around 5 MeV which lie well 
inside the spinodal line. 

As discussed in section 3 and 4, the spinodal regime 
does not deviate from the one obtained in the adiabatic 
limit. Although the growth rates for the spinodal insta- 
bilities are considerably reduced when realistic two-body 
collisions are taken into account, the decay times are still 
short enough such that spinodal decomposition occurs with- 
in the times the droplet remains around the turning point 
inside the spinodal line. Furthermore, while the density 
fluctuations grow, the restoring force of the radial flow 
is reduced, and hence the expansion tends to continue 
such that the spinodal decomposition is likely to complete. 
Therefore, we expect the qualitative features of multifrag- 
mentation as discussed in section 1.5 to survive for realistic 
dissipation. 



During the clustering process thermalization within 
the decaying droplet is so fast (relaxation time r rj 20 fm/c 
for T k 5 MeV) that statistical equilibrium is attained 
throughout the formation of fragments. This explains why 
statistical models are so successful in describing various 
aspects of multifragmentation. This may change as soon 
as the expansion is not completely slowed down and large 
flow velocities tear the droplet dynamically apart without 
allowing complete equilibration ||,|l(| . 

Recently, evidence of bulk fragmentation and spinodal 
decomposition have been reported in mass distributions 
for the collisions 155 Gd(36 MeV/u)+U and 129 Xe(32 MeV 
/u)+Sn, cf. Jnjjl^ . However, since mass distributions are 
dominated by statistical decay probabilities, the identifi- 
cation of spinodal decomposition turns out to be difficult. 
The situation may improve when the onset of spinodal 
decomposition is studied by looking into excitation func- 
tions. 
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